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We show that Bloch oscillations of ultracold fermionic atoms in the periodic potential of an optical 
lattice can be used for a sensitive measurement of forces at the micrometer length scale, e.g. in the 
vicinity of dielectric surface. In particular, the proposed approach allows to perform a local and 
direct measurement of the Casimir- Polder force which is, for realistic experimental parameters, as 
large as 10 -4 gravity. 
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The study of forces at small length scales consti- 
tutes one of the challenges of current frontier research in 
physics. The interest spans from fundamental issues, like 
the search of possible deviations from newtonian gravity, 
to the phenomenology of forces close to surfaces. Accu- 
rate investigations have been performed using a variety of 
methods, as for instance the techniques based on micro- 
cantilevers 0,0- The laser production and manipulation 
of micron-size ultracold atomic samples has provided a 
new tool for the investigation of forces at small length 
scales 0, 0- More recently, also trapped atomic Bose- 
Einstein condensates have been applied to the study of 
forces close to surfaces [H El • 

The possibilities offered by quantum degenerate gases 
have been enriched by the combination with periodic 
trapping potentials, as those produced by optical lattices. 
As early demonstrated H3 > this combination allows 
to study the general phenomenon of Bloch oscillations 
[Tl| . A key feature of this phenomenon is that the mea- 
surement of an oscillation frequency, which can be made 
very precise, can be turned into a direct measurement of 
a force. An important step in this direction has been the 
observation of temporally-resolved Bloch oscillations of 
a degenerate Fermi gas of atoms trapped in an optical 
lattice and subjected to gravity 01 ■ This has shown the 
feasibility of a high sensitivity force sensor with spatial 
resolution of the order of few microns. 

In this paper we explore the effect of a weak inhomoge- 
neous force close to a surface, such as the Casimir-Poldcr 
force ^| , on atomic Bloch oscillations (BO). We make a 
complete theoretical analysis for a Fermi gas trapped in 
an optical lattice in proximity of a dielectric surface, us- 
ing realistic parameters from an existing experiment |12| . 
We discover that this atomic Bloch oscillator realizes a 
powerful sensor for the detection of forces at small length 
scales. 

The physical system that we consider here consists of a 
sample of ultracold atoms trapped in a ID lattice aligned 
along gravity, as shown in Fig.^ The atoms are initially 
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FIG. f: Sketch of the physical situation considered in this 
paper: an ultracold sample of fermionic atoms is trapped in a 
vertical optical lattice in proximity of a dielectric surface and 
performs Bloch oscillations under the combination of gravity 
and of the Casimir-Polder force. 

cooled in a harmonic potential, and then transferred to 
the optical lattice. A harmonic horizontal confinement 
is assumed to be present, e.g. given by the same red- 
detuned lattice beams. As soon as the vertical harmonic 
confinement is switched off, the atoms start to perform 
BO in the lattice potential under the action of gravity. 
If now a surface is brought close to the atomic sample, 
additional forces between the atoms and the surface will 
affect the dynamics of BO. We will now investigate how 
such forces can be measured from the shift of the Bloch 
oscillation frequency as the distance between atoms and 
surface is varied. 

In particular, we consider here the Casimir-Poldcr 
force, which is exerted on an atom by the modified fluc- 
tuations of the electromagnetic field in the vicinity of a 
dielectric surface [l^. As a first step, we have worked out 
the simplified model in which the Casimir-Polder (CP) 
force Fcp is assumed to be spatially homogeneous. This 
is equivalent to assuming that the size of the sample along 
the vertical direction is small compared to the actual spa- 
tial variation of the force. In this case the quasimomen- 
tum q of the atoms simply evolves as 

hq = F = mg + F C p (I) 
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F being the total force acting on the atom. As q is de- 
fined only modulo the Bragg momentum qg = 2ir/\ (A 
is the wavelength of the optical lattice), each time q gets 
to the edge of the first Brillouin zone at qB, it reappears 
on the other side at — qs- The period of BO is thus 
Tb = iwh/FX. Following the evolution of BO for many 
periods one can infer their frequency, and therefore the 
Casimir-Polder force, with high accuracy. The fact that 
the frequency shift here is proportional to the force and 
not to its spatial gradient is an important difference of 
the present set-up with respect to previous work , and 
opens new possibilities for the measurement e.g. of ra- 
diation pressure forces close to a heated surface Q. In 
the effect of the CP force was investigated by look- 
ing at the shifts of the collective oscillation frequencies of 
a Bose-Einstein condensate confined in a harmonic trap. 
In this case the shift is proportional to the gradient of the 
force. Experiments in this direction are being performed 
by the group of E. Cornell at JILA |15j . 
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FIG. 2: Momentum distribution of a Fermi gas performing 
Bloch oscillations at different times within an oscillations pe- 
riod: solid lines are experimental data, dashed lines are the- 
oretical calculations including adiabatic switch-off of the lat- 
tice. E F /2S = 0.9. 

To observe BO one can study the momentum distri- 
bution of the atomic sample after a ballistic expansion. 
The typical evolution of the momentum is the one de- 
picted in Fig-121 At the initial time t — the momentum 
distribution is centered at k = 0, with symmetric lateral 
peaks at multiples of 2qs (lij- As time goes on, the comb 
of peaks moves rigidly with constant velocity, while their 
relative intensity follows the Fourier transform of a Wan- 
nier state [lTj as an envelope. In particular, when q(t) 
reaches the edge of the first Brillouin zone at t = Tg/2, 
the distribution contains a pair of symmetric peaks at 
±qs- At t = Tb, the distribution has regained the ini- 
tial shape. Note that the dynamics of BO in momentum 
space, being governed by the evolution equation JQ| for 
the quasimomentum, is independent of the lattice height, 
while the corresponding oscillations in real space have an 
amplitude of order S/F, where 28 is the width of the 
first energy band of the lattice. The presence of a strong 
known force, such as gravity, is therefore necessary to re- 
duce the amplitude of the oscillation and consequently to 
achieve a high spatial resolution in the measurement of a 
weak force. The use of a vertical geometry is helpful also 
to prevent the atoms from hitting the surface during the 



ballistic expansion. 

This discussion, originally developed for a single parti- 
cle, is immediately extended to the case of a many-atom 
gas simply by convolving the single particle prediction 
by the initial quasi-momentum distribution of the gas: a 
fundamental requirement to clearly follow BO is that the 
quasi-momentum spread of the sample must not be much 
larger than qs- This implies that the atomic tempera- 
ture must be at most of the order of 28. As the system 
size is limited by the spatial resolution we intend to ob- 
tain, a high density sample has to be used to increase 
the number of atoms and thus improve the signal. This 
motivates the use of a quantum degenerate gas. 

Although a Bose-Einstein condensate would certainly 
offer a much narrower initial quasi-momentum spread 
than a Fermi gas, it has been shown in recent experi- 
ments |18| to suffer from instabilities due to interactions, 
which quickly broaden the quasi-momentum distribution 
and make the observation of BO practically impossible 
in condensates at high density. As a matter of fact, for 
the rest of the paper we will mainly consider a Fermi gas: 
thanks to the absence of interactions this appears to be 
an excellent candidate for this kind of experiments |12| . 
For simplicity we will consider a Fermi gas at T = 0, 
whose quasi-momentum distribution is: 
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where E q is the energy of the Bloch state of quasi- 
momentum q, Ep is the Fermi energy, and uj r is the 
frequency of the radial harmonic trapping. As long as 
Ep is smaller than 28, the quasi-momentum distribution 
is localized in a region around q — and is zero outside, 
while the contrast slowly deteriorates for larger values of 
Ep. Notice that this behavior is different from the case 
of a purely one-dimensional system, where the contrast 
would be strictly zero as soon as Ep > 28 |l9j. 

The contrast can be improved by adiabatically switch- 
ing off the lattice potential at the end of the BO. In this 
case, the population of the state of quasimomentum q 
is projected onto the first Brillouin zone without being 
weighted by the Fourier transform of the Wannier state 
and at most two peaks are observed. In Fig.|3we report 
the result of a theoretical analysis taking into account 
the adiabatic switch-off of the lattice, which reproduces 
with excellent accuracy the experimental data taken as 
described in [l2|. 

Let us now study the effect of the CP force between the 
atoms and the dielectric surface. A general discussion of 
its properties has been recently given by , and we refer 
to that paper for the details of the calculations. In the 
simulations for the present paper, the potential plotted 
in Fig. 3 of has been used, which takes into account 
all the relevant regimes (van der Waals-London, Casimir- 
Polder, and thermal), in particular the thermal one 
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FIG. 3: Frequency shift of Bloch oscillations due to Casimir- 
Polder corce. The solid line is the prediction of Eq^Jbased on 
the thermal potential circles are the result of a complete 
numerical simulation as described in the text. 



which dominates at distances larger than the thermal 
wavelength hc/ksT of the photon. Here «o is the static 
atomic polarizability, e is the static dielectric constant 
of the material composing the surface, which is placed 
at z — 0, the z axis being oriented downwards. At the 
lowest level, one can approximate the spatially varying 
force with a homogeneous one equal in magnitude to the 
value of thermal force © evaluated at the center of the 
cloud. In this case, the BO period for an atomic cloud 
at a distance D from the surface is easily evaluated from 
cq[I]and uas a relative shift 



AT B _ F C p{D) _ 0.1728 
~T^~ ~ mg ~ D 4 



(4) 



The numerical value refers to the specific case of potas- 
sium atoms and a sapphire surface with eo = 9-4 and 
T = 300 K. In Fig.[3]we show the relative BO period shift 
as a function of the distance D, computed according to 
the approximated formula (J2J . Note how the effect of the 
Casimir potential ranges from 10 -4 to 10 -5 at distances 
D from 5 to 10 microns, which are realistic for a pos- 
sible experiment. This situation is encouraging because 
such level of sensitivity has been already demonstrated 
in a proof-of-principle experiment [l2T |. and can be fur- 
ther improved by a few orders of magnitude, as we will 
discuss later. 

It is then worth performing a deeper analysis tak- 
ing fully into account the spatial inhomogeneity of the 
Casimir force over the atomic cloud, its full z dependence 
beyond the large distance approximation (3) as well as 
the effect of a real-space motion of the cloud during BO. 
We have indeed performed an exact numerical calcula- 
tion for the non-interacting many fermion problem, in 
the combined potential of the trap, the gravity, the lat- 
tice and the complete CP force as predicted in Q. This 
has been done by first obtaining the single-particle or- 
bitals in the initial potential of the trap, the lattice, and 
the CP force, by then populating them according to a 
Fermi distribution, and by finally propagating them via 
the Schrodinger equation in the combined potential of 
the lattice, the gravity and the CP force. This proce- 



dure provides the density profile of the g clS, (IS well as 
its momentum and quasi-momcntum distributions at all 
times. From the center-of-mass of the quasi-momentum 
distribution, it is then possible to extract a theoretical 
prediction for the BO frequency shift. 




FIG. 4: Upper panel (a): momentum distribution of a Fermi 
gas at a distance D — 10 /im from the surface at beginning 
(dashed) and after 1000 periods of Bloch oscillations (solid). 
Lower panel (b): momentum distribution of a BEC after 2 
periods of Bloch oscillations. 

The results of such a calculation are shown in Fig.|3 
for realistic parameters: a situation has been considered 
in which 2 • 10 4 40 K atoms are prepared in a lattice 
with A = 873 nm, with a vertical extension of about 
4 /im, i.e. 8 lattice sites. To achieve a favourable ratio 
Ep/25 = 0.9 close to one we have chosen a lattice height 
s = 5, and a radial frequency of 10 Hz. The vertical ex- 
tension in the lattice is determined by the confinement 
of the initial harmonic potential of a trap, which for this 
specific case has a frequency of 400 Hz. Justification for 
the slight difference with respect to the analytical for- 
mula (0J can be found in several effects that have been 
neglected in the previous discussion, in particular the 
spatial displacement (on the order of 1 /Ltm in the direc- 
tion of going further from the surface) of the atomic cloud 
while performing BO, the difference between the approx- 
imate thermal potential © and the exact one of and 
finally the spatial inhomogeneity of the CP force across 
the atomic cloud. The inhomogeneity of the force can 
have a further consequence: the frequency shift of BO is 
in fact given by an average of the force over the cloud 
profile, and a sort of inhomogeneous broadening is to be 
expected. As this latter effect may affect the accuracy 
of the measurement, in particular after the large num- 
ber of oscillations which are required for a high-precision 
measurement, we have investigated the broadening of the 
momentum distribution after a large number of BO pe- 
riods. The results, reported in Fig.0Ji, show that the 
broadening is not yet significant after 1000 periods. 

For the sake of completeness, a complete analysis has 
been also performed for the bosonic case by numerically 
solving the ID Gross-Pitaevskii equation with trap and 
atom number parameters analogous to the Fermi case 
(that is with a peak density of n = 4 ■ 10 13 cm -3 ). The 
results are in agreement with the experimental obser- 
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vations of 12]: the quasi-momentum distribution dra- 
matically broadens as soon as the quasi-momentum of 
the Bose condensate enters the dynamically instable re- 
gion [2(j and after only a few BO periods, it is not pos- 
sible to follow the oscillations any longer (Fig.|3J>). To 
overcome this problem, one could think of setting the 
scattering length a to zero by means of a Feshbach res- 
onance [21|. As the characteristic time for the onset of 
the instability is of the order of r, ns t « h/gn max |20| . 
9 n max = 47r?i 2 an max /m being the mean-field interac- 
tion energy, an upper bound on the absolute value of a 
is set by the fact that Ti ns t has to be much longer than 
the actual duration of the measurement. For instance, 
a measurement lasting 100 BO periods, performed with 
the atomic density here considered, would require a re- 
duction of the scattering length of a factor of the order 
of 10 3 . 

Let us now discuss the sensitivity that can be real- 
istically achieved in an experiment. In |l2^ we demon- 
strated an overall sensitivity ATb/Tb = 10~ 4 , by follow- 
ing about 100 BO periods of 40 K atoms in a lattice with 
A = 873 nm. This sensitivity was limited by various 
sources of decoherence due to the imperfect nature of 
the optical lattice. Examples are the phase noise due 
to mechanical vibrations of the retroreflecting mirror, or 
longitudinal inhomogeneous forces due to the focusing of 
the lattice beams. Most of these decoherence sources can 
be reduced in an optimized experiment, thus extending 
the number of oscillations observable with high contrast. 
The typical lifetime of the atomic sample of a few sec- 
onds will however limit the realistic number of observable 
BO periods to the order of 10 3 . Since the sensitivity on 



the measurement of a single BO period can be pushed to 
10~ 3 -^10 -4 , as typically achieved in this kind of experi- 
ments |22j , one can expect to obtain an overall sensitivity 
AT B /T B =10- 6 ^10- 7 . 

This sensitivity would certainly allow to study a weak 
and inhomogeneous force like the Casimir-Polder force 
already at distances larger than 5 /xm from a dielectric 
surface. This level of accuracy might also open to the 
investigation of possible deviations from Newton's grav- 
itational law at short distance [2^]. In order to improve 
the probing of the force at shorter distances one should 
reduce the width of the sample by working, for example, 
with a stronger vertical confinement and with a shorter 
lattice periodicity. The latter would indeed work in favor 
of a higher resolution, without affecting the sensitivity 
because, as already discussed, the measurable quantity q 
depends only on the applied force. 

In conclusion we have discussed how Bloch oscillations 
can provide a sensitive tool for the study of forces at 
the micrometer scale, e.g. the Casimir-Polder force be- 
tween atoms and dielectric surfaces. The frequency shift 
of Bloch oscillations is proportional to the applied force 
and at a distance of the order of a few microns turns out 
to be well within the sensitivity of realistic experiments. 
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